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LIPSCHITZ-KILLING CURVATURES OF SELF-SIMILAR 
RANDOM FRACTALS 

M. ZAHLE 

Dedicated to Herbert Federer 

Abstract. For a large class of self-similar random sets F in R d geometric 
parameters Ck(F), k = 0, . . . , d, are introduced. They arise as a.s. (average or 
essential) limits of the volume Cd(F(e)), the surface area Cd—i{F(e)) and the 
integrals of general mean curvatures over the unit normal bundles Ck(F(e)) 
of the parallel sets F(e) of distance e rescaled by e D ~ k as e — > 0. Here D 
equals the a.s. Hausdorff dimension of F. The corresponding results for the 
expectations are also proved. 



0. Introduction 

Self-similar sets in Euclidean space R d in the sense of Hutchinson [11] have 
intensively been studied in fractal analysis and geometry. Their probabilistic coun- 
terparts, the so-called (stochastically) self-similar random sets, were introduced 
independently by Falconer [2], Mauldin and Williams [17] (1986), and Graf [8] 
(1987). Concerning the geometry of these random sets the literature up to now 
deals mainly with exact Hausdorff dimensions, multifractal spectra and associated 
measures. 

In the present paper we will introduce a system of d + 1 (random) geometric pa- 
rameters which allow to distinguish between self-similar random sets F with equal 
Hausdorff dimension D, but different geometric and topological features. This 
continues the work of Winter [25] who first investigated such parameters for deter- 
ministic self-similar sets with the open set condition and polyconvex neighborhoods. 
Here we use the construction model from [17] for such random F satisfying the open 
set condition for a fixed deterministic open set with closure J. Our parameters are 
in close relationship to Federer's curvature measures for sets of positive reach [3]. 
For this we assume additionally that with probability 1 for Lebesgue almost all 
r > the parallel set F(r) of amount r has a Lipschitz boundary and the closure of 
its complement F(r) has positive reach. (We conjecture that this is already guar- 
anteed by the strong open set condition on int(J). At least this can be checked for 
many examples which do not possess polyconvex parallel sets. For d < 3 it follows 
from a general result of Fu [6]). In this case the random sets F(r) are Lipschitz 
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d-manifolds of bounded curvature in the sense of [22] and their Lipschitz-Killing 
curvatures are determined by those of Federer: 

C k {F{r)) = {-l) d - 1 - k C k (W)) ,fc = 0,...,d-2. 

They are completed by the volume Cd{F{r)) and the surface area Cd-i(F(rj) with- 
out the additional assumption. The equalities for the total curvatures are localized 
to (signed) curvature measures. Our main result (see Theorem 2.3.8 and Corol- 
lary 2.2.4) is the following: Under an additional assumption on the model (only if 
k < d - 2) the limit 

l 

8 

and the almost sure limit 

l 

^jEff J >" fc <W(s)) ^ 

exist. Moreover, the first limit has an explicit integral representation and the sec- 
ond one is a random multiple of the first one. The multiplier is the inverse of a 
martingale limit related to the contraction ratios of the generating system of ran- 
dom similarities, i.e. it does not depend on k. If the logarithmic contraction ratios 
are non-lattice in a sense, then the average limits may be replaced by the ordi- 
nary limits as e — > 0. In general, the limits define random fractal Lipschitz-Killing 
curvatures of F, if they do not equal zero. For vanishing limits the correctness of 
the scaling exponents D — k has to be checked in order to give such a curvature 
interpretation. 

For the special case k = d this concerns the (average) Minkowski content, and 
the results were proved by Gatzouras (see [7] and the references therein to related 
work). We adopt his idea to apply the classical renewal theorem for the expecta- 
tions and a renewal theorem of Ncrman [18] for branching random walks in order 
to derive the almost sure limits. 

At the same time we get an extension of Winter's results for the total curvatures in 
the deterministic case. However, Winter also proved such limit relationships for the 
corresponding curvature measures. The fractal versions in his case are all constant 
multiples of D-dimensional Hausdorff measure on F. It will be shown later that 
this remains valid under our conditions. In order to introduce random fractal cur- 
vature measures additional work is needed, in particular, related to exact Hausdorff 
measures of such random sets in the sense of Graf, Mauldin and Williams [9] . 
The paper is organized as follows. 

In Section 1 we summarize some background from classical singular curvature the- 
ory. The Appendix provides an auxiliary new result result in this direction - an 
estimate for the Lipschitz-Killing curvatures of sufficiently large parallel set of an 
arbitrary compact subset of M. d , which is related to well-known isodiametric in- 
equalities in the convex case. 

Section 2.1 provides the random iterated function scheme and the limit set F by 
means of random trees and an associated branching random walk. 
In Section 2.2 we follow Gatzouras [7] and present a slight extension of the corre- 
sponding renewal theorem for branching random walks. 



LIPSCHITZ-KILLING CURVATURES OF SELF-SIMILAR RANDOM FRACTALS 3 

Its application to fractal curvatures is prepared in Section 2.3, which contains the 
formulation of the main results. 

The key proofs are given in Section 3 using the above mentioned theorem from the 
Appendix. 

1. Background from classical singular curvature theory - curvature 

measures of parallel sets 

The geometry of classical geometric sets in R d may be described by certain pa- 
rameters forming complete systems of Euclidean invariants in the following sense: 
It is well known from convex geometry (Hadwiger's characterization theorem) that 
every set-additive, continuous and motion invariant functional on the space of con- 
vex bodies is a linear combination of Minkowski's quermass integrals. The latter 
are also called intrinsic volumes. For smooth boundaries they agree - except the 
volume - with the integrals of the elementary symmetric functions of principal cur- 
vatures, the so-called integrals of (higher order) mean curvatures of smooth man- 
ifolds. Federer [3] unified and extended the intrinsic volumes and the integrals of 
mean curvatures with tools of geometric measure theory: He introduced curvature 
measures of sets with positive reach by means of a Steiner polynomial for the vol- 
ume of parallel set. (This is related to Weyl's tube formula in differential geometry, 
where the mean curvatures arise as traces of certain powers of the Riemannian 
curvature tensor.) Explicit representations of Federer's curvature measures in form 
of integrating Lipschitz-Killing curvature forms or symmetric functions of gener- 
alized principal curvatures over the associated unit normal bundle were given in 
[28]. Starting from the 1980s (and earlier in convex geometry) up to now additive 
extensions and other versions of these curvature measures for classes of singular sets 
in M. d have been studied with methods of geometric measure theory and algebraic 
geometry. Nowadays the notion Lipschitz-Killing curvature measures is used. 
Localizations of Hadwiger's characterization theorem to curvature measures of con- 
vex sets where given by Schneider [23]. In [29] this was used together with an 
appropriate notion of continuity as approximation tool in order to generalize these 
characterizations to large classes of singular sets. At the same time this served as 
a motivation to study such curvature properties also for fractals. 
The first results in this direction were obtained by Winter [24], [25], who worked 
out essential tools for investigating self-similar fractal sets in K d under this point of 
view. Under the usual open set condition and the additional assumption of poly- 
convex neighborhoods he solved the corresponding problems completely. In view of 
the Gauss-Bonnet theorem the notion of fractal Euler number investigated before 
in Llorente and Winter [14] under some additional assumption may be considered 
as a special case. 

In order to avoid the condition of polyconvex parallel sets and to extend such cur- 
vatures to random fractals we suggest now another approach. For certain classes 
of compact sets K C M. d (including many classical singular sets) it turns out that 
for Lebesgue-almost all distances r > the parallel set 

(1) K{r) := {x e R d : d(x, K) = min \x - y\ < r} 

possesses the property that the closure of its complement 



(2) 



K(r) := K{r) c 
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is a set of positive reach with Lipschitz boundary. A sufficient condition is that 
r is a regular value of the Euclidean distance function to K (sec Fu [6, Theorem 

4.1] together with [21, Proposition 3]). In this case both the sets K(r) and K(r) 
are Lipschitz d-manifolds of bounded curvature in the sense of [22], i.e., their k-th 
Lipschitz- Killing curvature measures, k = 0,1,..., d— 1, are determined in this 
general context and agree with the classical versions in the special cases. Moreover, 
they satisfy 

(3) C k (K(r),-) = (-l) d - 1 - k C k (K(r),.). 

Therefore the C k (K(r), •) are signed measures with finite variation measures 
C k mr {K{r) 1 •) and the explicit integral representations are reduced to [28] (cf. [22, 
Theorem 3] for the general case). In the present paper the normal cycle represen- 
tation is needed only in the Appendix. We will briefly mention the normal bundle 
and the current construction at the end of the section. (The reader not familiar 
with the corresponding geometric integration theory may skip the current theoret- 
ical part.) Here we will list the main properties of the curvature measures for the 
parallel sets as above which will be used repeatedly: 

Cd-i(K(r), •) agrees with (d — 1)- dimensional Hausdorff measure on the 

boundary dK(r). The latter is a bounded measure for all r > and all com- 
pact sets K to which it is applied below. Therefore we use this notation in any 
case. Furthermore, for completeness we define Cd(K(r),-) as Lebesgue measure 
restricted to K(r). The total measures (curvatures) of K(r) are denoted by 

(4) C k (K(r)) :=C k (K(r),R d ), k = 0,...,d. 

By an associated Gauss-Bonnet theorem (see [21, Theorems 2,3]) the Gauss curva- 
ture Co{K(r j) coincides with the Euler-Poincare characteristic x(A(r)). 
The curvature measures are motion invariant, i.e., 

(5) Ck(g(K(r)),g(-)) = C k (K(r), •) for any Euclidean motion g , 
they are homogeneous of degree k, i.e., 

(6) C k (XK(r),X(-)) = X k C k (K(r),-), A >0 , 
and locally determined, i.e., 

(7) C k (K(r), (•) n G) = C k (K'(r'), (•) n G) 

for any open set G C R d such that K(r) n G = K'(r') n G, where K(r) and A'(r') 
are both parallel sets where the closures of the complements have positive reach. 

We now summarize some facts about sets with positive reach needed below: 
Recall that reach A of a set A C Mr is defined as the supremum over all s > 
such that every point x in the s-parallel set of A there is a unique point Ilxa; G A 
nearest to x. The mapping II^ (on its domain) is called the metric projection onto 
A. For a set X of positive reach the unit normal bundle is defined as 

nor A := {(x, n) E R d x S^ 1 : x E A, n e Nor(A, x)} 

where Nor(A, x) is the dual cone to the (convex) tangent cone of A at x. 

If additionally nor A n p(nor A) = for the normal reflection p (x,n) [x, —n), 

then A is full dimensional with Lipschitz boundary (see [21, Proposition 3]). 
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For general X with reach X > there is an associated rectifiable current called the 
unit normal cycle of X which is given by 

N x (ip):= J (a x {x,n),^(x,n))n d - 1 {d{x,n)) 

nor X 

for an appropriate unit simple (d — l)-vector field ax = &i A . . . A a^-i associated 
a.e. with the tangent spaces of nor and for integrable differential (d — l)-forms 
(p. In these terms for k < d — 1 the curvature measure may be represented by 

C k {X,B) =iV X Ll BxR£i ((p fe ) = J (ax(x,n),ip k {n))U d -\d(x,n)) 

norXn(SxR d ) 

for any bounded Borel set B C R d , where the fc-th Lipschitz-Killing curvature form 
tfk does not depend on the points x and is defined by its action on a simple (d— 1)- 
vector r] = r]i A . . . A f]d-i as follows: Let n (y, z) := y and iri(y,z) := z be the 
coordinate projections in R d x R d , e[, . . . , e' d be the dual basis of the standard basis 
in R d and Ok the surface area of the fc-dimensional unit sphere. Then we have 

(V,<Pk(n)) := Ol\ (7r £l ?7i A...Avr ed _ 1 77 (i _iAn,e / 1 A...Ae^). 

£.e{o,i},E£.=ii-i-i! 

2. Self-similar random fractals - the model and statement of the 

main results 

2.1. Random recursive constructions and associated branching random 
walks. We briefly describe the random recursive construction model introduced in 
Mauldin and Williams [17] (and independently in Falconer [2] and Graf [8] with 
different methods). Additionally we use some ideas from Gatzouras [7] for rela- 
tionships of associated random functions to branching random walks. Let Sim be 
the set of contracting similarities and J a nonempty compact subset of of R d with 
J = int(J). Our basic object is a random element S :— {Si, . . . , S„} with Si € Sim 
if v > 0, and S := {id} if v = 0, where v is a random variable with values in 
No := {0} U N. We suppose that S satisfies the open set condition (briefly (OSC)) 
with respect to int(J): 

V 

(8) |J 5i(int(J)) C int(J) and Sj(mt(J)) n Sj{mt(J)) = 

i=l 

with probability 1. Denote the corresponding probability space by [f2°, P ]. S 
is also called random iterated function system, briefly RIFS. The corresponding 
random fractal set is introduced by means of the code space E := N N and a random 
Galton-Watson tree in the set of all finite sequences := {0} U U^Li 
For a = <7i . . . ak , t = t\ . . . ti £ we write \a\ := k for the length of a, u\i :— 
o~i . . . o~i, ijk, for the restriction to the first i components, and or := o~\ . . . o~kT\ ■ ■ ■ V, 
for the concatenation of a and r. (We will use analogous notations for infinite a e S, 
resp. t e S.) By convention, Oct = a. 

For each a £ £* let [f2 CT , So-, P<r] be a copy of the above probability space. The basic 
probability space for the random construction model is the product space 

(9) [fi,ff,P] := (g)[0 CT ,^,P ff ]. 
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On this space a family of independent identically distributed RIFS 
{S a = {Sal, ■ ■ ■ , So-ivHo-eE, (where S a = {id} if v„ = 0) with i.i.d. random 
numbers {i/ a } aeS , is then determined. Let $ n be the cr-algebra generated by all 
S a and v a with \a\ < n. For brevity we write 





= Sa\l O S a \2 ° ■ ■ ■ ° S a \\a\ , 


r a 






= Up(Sa) = ra\ir a \ 2 ■ ■ ■ r a \\a 


Ja 


= S a J, 


K a 


= KnJ a , 


K a 


= S-\K)nJ, 


\K\ 


= diam(isr) , 



a E S*, for any compact set K. Note that K a is a random compact subset of J in 
the sense of stochastic geometry (measurable with respect to the Borel cr-algebra 
given by the Hausdorff distance). 

Set 7o := {0} and define inductively T n +\ '■= 0, if T n = 0, and 
Tn+i := {cri: a € T n , v a °) 1 < * < v <y} , 

if T n ^ 0- Then 

OO 

T:= [jT n 

is the population tree of a random Galton- Watson process. T n represents the family 
of individuals in the n-th generation with ancestor 0. The boundary of T is defined 

by 

dT := {a £ S : a\n £ T, n £ N} . 
In the sequel we consider the supercritical case with 

(10) l<Ei/<oo, 

where this boundary is nonempty with positive probability (see the classical liter- 
ature on branching processes). The random compact set 

OO 

(11) F:= p| (J J a 

n=l aeT n 

is the associated self-similar random set. F is the image of the boundary dT under 
the random projection 

tt : a M> lim S a \ n (x ) 

n— foo 

for an arbitrary starting point xq £ M. d . By construction the random set F is non- 
empty with positive probability. Its stochastic self- similarity property reads here as 
follows: 

(12) F= |J StF*) 

StGS 

where the random sets F l , i £ No, are independent, have the same distribution as 
F, and the random clement S = {Si, . . . , S„} with contracting similarities Si, if 
v > 0, and S — {id}, if v = 0, is as above and independent of the F l . 
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It is well-known that with probability 1 (briefly w.p.l) the Hausdorff dimension D 
of the self-similar random set F is uniquely determined by the equation 

(13) E(5>f)=l. 

Under the additional assumption that P(F fl int( J) ^ 0) > 0, the so-called strong 
open set condition for the open set int(J) from OSC (8), this has been proved in 
[17], [2] and [8]. Following Patzschke [19] it remains true supposing only OSC. 
We are interested in curvature properties of the random fractal set F . The main 
tool will be approximation by parallel neighborhoods of small distances using the 
Lipschitz-Killing curvature measures from Section 1 and suitable rescalings. (For 
the case of deterministic similarities (Si, . . . , Sn) and the assumption of poly convex 
neighborhoods for the deterministic self-similar set F the corresponding notions and 
results have been worked out in the Thesis of Winter, see [25].) The related problem 
for the Minkowski content solved in [7] may be considered as a marginal case. An 
important tool will be again the branching random walk defined by the 

recursive formula 

(14) W <r = W <r | (W _i)+Inr- 1 , 

if a e T, and W a := oo, if er e E*\T. (Recall that r a denotes the contraction ratio 
of the similarity S a for a e T, and by convention ro = 1.) In particular, Wq = 0. 
As in [7] the classical renewal theorem essentially used for the above problems in 
the deterministic case will be replaced by an associated stochastic version. 

2.2. Renewal theorem for branching random walks. We refer to Gatzouras 
[7, Section 3.2]. First consider the associated nonnegative martingale 

(15) M n := e ~ DW °= E f ?< »^ ' 

with respect to the filtration {3n}n>o- (Note that by (13) EMi = 1.) According 
to the martingale convergence theorem the limit 

(16) Moo := lim M n 

n— >oo 

exists w.p.l. The next theorem is due to Biggins [1]; see Lyons [15] for a conceptual 
proof using Lyons, Pemantle and Peres [16]. 

Theorem 2.2.1 (Biggins). The following are equivalent: 

(i) E(Miln+Mi) < oo 

(ii) ¥{Moo - 0) < 1 

(iii) P(Mco > | non- extinction) = 1 
(vi) EMqo = 1 

(v) M„ Moo in L 1 . 

Recall that the underlying probability space [0,#, P] = [fi<r , da , Per] is a 

product space. For each reE, define the shift operator 8 T : Q — > Q by 

(17) (OtU)o- ■= LO Ta . 

Besides the branching random walk {W a } ae j- from (14) we now consider a basic 
stochastic process Y satisfying the following. 
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Condition 2.2.2. Y = {Yt : t G R} is a real-valued measurable stochastic process 
on [£1,3, P] vanishing for t < 0, which is continuous a.e. with probability \, i.e., 
there exists a set C £ 93 (R) $ such that J Q J R lc {t, w) dt P(duj) = and for any 
(t,iv) G C the function Y^(uj) is continuous at point t. 

The process Y induces a family of i.i.d. copies defined by 

Y t a (u):=Y t (O a w), <J6S t . 

Then we can introduce the branching process associated with W and Y in the sense 
of Jagers [12]: 

(18) Z t := Y t-w ■ 

<t£T 

We are interested in the limit behavior of the process e~ Dt Z t as t — > oo. Here the 
measure 



!(■)(! lnr * 



(19) /i-Ef^l^Wi)) =e(E : 

on R plays a crucial role. Note that n(M.) = and recall the assumption 1 < Ev < 
oo. Denote 

(20) \(D):=E(j2 W * e ~ DWi ) =E(2l lnr 'l^ 

\ieTi / VieTi 



The following renewal theorem for the branching process Z is essential for our pur- 
poses. (Recall the random martingale limit Mqo from (16)-) Note that the essential 
limit as t — > is meant w.r.t. Lebesgue measure and is defined as the common value 
of limesssup and limessinf, if the latter coincide. 



Theorem 2.2.3. Suppose that the process Y satisfies Condition 2.2.2 and there 
exists a non-increasing integrable function h : [0, oo) — > (0,oo), such that 

(21) E ess sup - ' 1 <oo. 

V t>o h{t) J 

(i) [Nerman] If the measure \i is non-lattice then 

oo 

esslime-^E^) = f e - Ds E(Y s )ds 



and 

oo 

ess\ime- Dt Z t = [ e- Ds E(Y s )ds w.p.l . 

t^oo J 


(ii) [Gatzouras] If the measure /i is lattice with lattice constant c, then for 
Lebesgue-a.e. s <G [0, c) we have 

oo 

lim e- D{s+n ^E{Z s+nc ) = — Ve- fl ( 8+ " c 'E(i; +BC ). 

\(D) ^ 
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and 

1\/T 00 

lim e- D ^Z s+nc = ]T e- Dis+n MYs +nc ) w.p.l . 

\(D) ^ 

Assertion (i) is shown in Nerman [18] for the case of Skorohod-regular processes 
(which are a.e. continuous w.p.l). The lattice case (ii) is derived in [7] from Ner- 
man's proof of (i). Note that these proofs remain valid under our Condition 2.2.2 
and (21) when considering essential limits. The convergence of the expectations 
may be considered as a special case of the classical renewal theorem for deter- 
ministic functions. (Feller's proof in [5] works also for the essential limits.) A 
straightforward consequence of (ii) is the following. 



Corollary 2.2.4. Suppose that the conditions of Theorem 2.2.3 (ii) are satisfied 
and 

c 

(22) J sup ( e - D ( a+nc >E \Z s+nc fj ds < 00 . 



Then we have 

T 00 



and 

T 00 


Proof. We show the almost sure convergence. The arguments for the expectations 
are similar. 

The branching process Z is measurable with respect to 03 (R) (g> $ and therefore 

£?:=((«, w)G [0,oo) xO: lim e - D (' +nc >Z a+nc (w) - L s ) e 03(R) O 3, 
where L s denotes the right hand side of (ii). Theorem 2.2.3 (ii) and Fubini yield 



= J j l E c(s,Lu)F(dLo)ds = j j l E c(s,uj)ds¥(dw), 

on no 

i.e., w.p.l we get lim n ^ oc e~ D( ~ s+nc ^Z s+nc = L s for a.e. s e [0,c). Taking into 
regard (22) (where we could omit the expectation for convergence w.p.l) and dom- 
inated convergence we infer 



c c 

lim / e- D ^ s+n ^Z s+nc ds= f L s 



ds 





w.p.l. Instead of the sequence on the left hand side we may also take its arithmetic 
means which converge to the same limit. The latter is equal to the right hand side 
of the assertion as an easy calculation shows. Similarly, the limit of the arithmetic 
means is the same as the limit of the Cesaro means on the left hand side of the 
assertion. □ 
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2.3. Application to fractal curvatures. We now turn back to the self-similar 
random set F from (11). Let JC be the space of non-empty compact subsets of 
our primary compact set J. Q3(/C) denotes the Borel d-algebra with respect to the 
Hausdorff distance djjoiiK. Wc further consider the space F of closed subsets of R d 
provided with the hit-and-miss topology (generated by the sets {A e T : A(lO ^ 0} 
and {A e T : A n C = 0} for open O and closed C) and the Borel-tr-algebra 23 (J 7 ). 
The topology restricted to JC is generated by the metric dn ■ 

(Recall the notations (1) and (2) for parallel sets and the closures of their comple- 
ments. It is easy to see, that the mapping 

(r, K) K(r) from [0, oo) x JC to T 
is (<8([0,oo)) <g> Q5(/C) , s B(J"))-mcasurablc. 

Let VIZ be the space of subsets of R d with positive reach (cf. the end of Section 1). 
According to [27, Proposition 1.1.1] we get VTZ € 25 (J 7 ). Moreover, the mapping 
X i y noiX from VIZ into the space of closed subsets of R d x R d is measurable. 
(Recall the definition of the unit normal bundle nor X from Section 1.) This implies 
the following measurability property (for the normal reflection p : (x, n) ^ (x, — n)): 

Lemma 2.3.1. 

Reg : = j (r, K) G [0, oo) x JC : K{r) e VIZ, nor K{r) n p(nor K(r)) = j 

€ «8([0, oo)) ® «8(/C) . 

The elements of i?eg will be called regular pairs. 
In Fu [6, Theorem 4.1] it is shown that in space dimensions d < 3 for any compact 
set K there exists a bounded exceptional set TV of Lebesgue measure such that for 
any r ^ N the set K (r) has positive reach and nor K (r) n p(nor K{r)) = 0. More- 
over, if r > y/d/ (2d + 2)\K\ both these assertions hold for any space dimension d. 
This basic result, the preceding lemma and Fubini imply the following property of 
the self-similar random set F: 

Corollary 2.3.2. In space dimensions d < 3 w.p.l for Lebesgue-a.e. r > the 
pair (r, F) is regular in the above sense. 

In higher dimensions we shall formulate this as a main geometric condition on 
the self-similar random set F: 

Definition 2.3.3. The self-similar random set F is called regular if w.p.l for 
Lebesgue-a.e. r > the pair (r, F) is regular. 

For regular pairs (r, F) the Lipschitz-Killing curvature measures of F(r) arc 
determined by means of those for F(r) (cf. Section 1): 

(23) C k (F(r),-) = (-i) d - 1 - k C k (Ff),-) 7 fc = 0,...,d-l. 

According to Fu's result mentioned above this holds, in particular, for any realiza- 
tion of F and r > r^ where 

r d := y/d/{2d+2)\J\. 
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For the exceptional pairs we set 

C k (F(r),-) :=0. 

Then the Ck(F(r), •) are random signed measures. The corresponding measurabil- 
ity properties follow from [27, Theorem 2.1.2] and Lemma 2.3.1: If the self-similar 
random set F is regular, the mapping (r, F) i-> C(F(r), B) for any Borel set B e K d 
is product measurable. Moreover, one obtains the following continuity property: 

Lemma 2.3.4. For any (tq,K) € i?e<? and r — > ro > we /lave (r, if) <G iieg and 
ffte measures (7fc(if (r), •) weakly converge to (7fc(if (ro), •). 

The proof follows from [21, Propositions 1 and 3] and [22, Proposition 6]. (Re- 
call that if (r, if) £ Reg both the sets if (r) and if (r) are Lipschitz d-manifolds of 
bounded curvature in the sense of [22].) 
This implies the following auxiliary result. 

Corollary 2.3.5. If the self-similar random set F is regular, then the random 
function $(r) := Ck{F(r)) is a.e. continuous w.p.l. 

Remark 2.3.6. If k = d this property remains valid for general compact sets and 
all r > 0, i.e. we need not restrict to the class VIZ. For k = d — 1 the set VTZ can 
be replaced by the set of d — l)-rcctifiable closed subsets of R d in the sense 

of Federer [4] . The above measurability and continuity properties for this case are 
treated, e.g., in [26]. (Here the weak convergence of Cd-i(K(r), •) to Cd-i(K(r ), •) 
as r — > r is a well-known result from geometric measure theory.) 

In the Appendix it will be shown that parallel sets of distances greater than 
\/2| J| are nice sets concerning their regularity properties. Therefore in the sequel 
we fix an arbitrary 

R>V2. 

It turns out that the relevant formulas below do not depend on the choice of R. 
In order to formulate the remaining conditions on F and to apply the renewal 
theorem we now turn back to the tree construction from Section 2.1: 
For fixed k G {0, 1, . . . d} we consider the basic stochastic processes 

(24) Z t := (R\J\)- k e kt C k (F (ii| J|e-*))l [0>oo) (t) 
and 

(25) Yt ~Z t -Y, Z l Wi ■ 

Iterating the last equation we get for Z t the branching process representation (18): 

Z t = Y t-w a ■ 

crGT 

(Recall the notation X^(uj) — X t (9 a (uj)), u e S,, for a process X and the shift 
operator (17) on the probability space (9).) Here we have in mind the scaling 
property and the additivity of the Lipschitz-Killing curvature measure C k in order 
to obtain the conditions on the process Y for the renewal theorem. The continuity 
condition 2.2.2 will follow from the regularity of the self-similar random set F if 
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k < d — 2 or from the rectifiability of the boundaries of its parallel sets if k = d — 1. 
For the boundedness condition (21) we will use an additional assumption if k < d—2, 
which is formulated in the language of random stoppings in the tree construction: 
For any r > we define the random subtree 

(26) T(r) := {a e T : R\J\ f a <r< R\J\ f CT | khl } . 
This is a so-called Markov stopping on our probability space. 

Recall the above notations F" = S~ 1 (F) n J = F o 6 a for the shift operator 9 from 
(17). Then we obtain from the representation F = U<xer(r) F a wrtn F a = F n J 
the stochastic self-similarity 

(27) F= |J Sa-(F a ) , 

<reT(r) 

where i ;icr , it G S*, are i.i.d. copies of F and the random similarities S„, ff £ E„ 
and the random subtree T(r) are as before and independent of the F° '. At the 
same time we get for any r > the representation 




For any a G T(r) we have > R\J\ and F a C J. Then Theorem 4.1 in the 
Appendix implies that 9(F cr (^-)) is a Lipschitz (d — l)-submanifold. Since the 
number of a G T(r) is finite, we have shown the following: 



Lemma 2.3.7. With probability 1 for any r > the boundary of the random set 
F(r) is {U d - 1 ,d-l)-rectifiable. 

(For the deterministic case see Rataj and Winter [20].) 
Furthermore, we introduce the subset of those words a in T(r) for which the set 
F a (r) = (F fl Ja)(r) has distance less than r to the boundary of the first iterate 
SJ := UieTi ^ i °^ ^ e basic set J under the random similarities: 

(28) T b (r) := {a G T(r) : F CT (r) n (5J) c (r) ^ 0} . 

We now can formulate the mam theorem of the paper. 
Recall the martingale limit = lim^oo ^2 ae7 - n from (16), the measure 

V : = E (E 4e r 1 !(■)(! lnr ^l)) from ( 19 )' and the constant \{D) = E (£igTi I In n\ if) 
from (20). 



Theorem 2.3.8. Let k G {0,1,..., d} and F be a self-similar random set with 
basic space J satisfying the following conditions: 
(i) 1 • .:',(' x . 

(ii) the strong open set condition, i.e., the open set condition (8) and 

¥(Fr\mt(J) ^ 0) > 0, 

(iii) F is regular in the sense of Definition 2.3.3, if k < d — 2, 

(iv) fork<d-2, 



E ess sup sup r- k Cr r {F(r),d(F rJ (r))nd( M F a , (r)) < oo . 

\0<r<fl|Jk€T 6 (r) V ^ 
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Set 

R k (r) :=C k (F(r))-J2Mo,R\J\r t ](r)C k (F t (r)), r >0 . 

i£Ti 

Then we have the following. 
(I) 

R\J\ 



ess 



Kms D - k E(C k (F(s))) = ^ J r D - k -'E(R k (r))dr 



and 

R\J\ 



C k {F):= G ssY T e D - k C k {F(e)) = ^ j r^-^E^r)) dr w.p.l. 

o 

provided the measure fj, is non-lattice. 



(II) 

n lim o e( fe - I3 )( s+ " c )E(c fc (F( e -( s +" c )))) = J e (k-D)(s+mc) E ^ fe ( £ -(.+™o))j 

' m=0 

/or a.e. s€ [0, c) and 

A/f °° 

Urn e ( fc -^)(-+™ c )c fe (F(e-( s +" c ))) = ^ e< fe -°X i+mc >E(i? fc (e-< s+mc >)) 

^ ' m=0 

for a.e. s G [0, c) w.p.l, 

provided the measure \x has lattice constant c. 

Remarks. For k = d, i.e. for the (average) Minkowski content in (I), this theorem 
is due to Gatzouras [7]. 

Recall that for d < 3 the regularity (iii) holds always true. For polyconvex neigh- 
borhoods this remains valid for general d. 

We conjecture that in the deterministic case the strong open set condition implies 
(iii) . Under the additional assumption of polyconvex neighborhoods (iv) is proved 
implicitly in the Thesis of Winter [24] , [25] . For many deterministic examples with 
non-polyconvex neighborhoods, e.g. the Koch curve, the above conditions can be 
checked using their local structure. 

For the general deterministic case and k = d — 1 the limits are derived in Rataj and 
Winter [20]. Moreover, these authors show for the case of non-arithmetic logarith- 
mic contraction ratios the equality 

C d - 1 (F) = (d-D)C d (F). 

This supports our conjecture that like in the classical smooth case the parameters 
C k (F) with k larger than the integer [D] + 1 do not provide additional geometric 
information. (In the classical case they are all multiples of the Minkowski content.) 

Corollary 2.3.9. Under the conditions of Theorem 2.3.8 we get the average limits 
l R\J\ 
iToUn^ I ^ D ^ E (C k (F(e))) e^de = ^-L_ J r D -^E(R k (r)) dr 
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and 

1 R\J\ 
l^>\^l\J £D ^ Ck{F{£))£ ^ d£= ^E) J rD ~ k ~ lE ( R k(r))dr w.p.l. 

8 

The notion of fractal curvatures. In the non-lattice case these average limits 
agree with the ordinary limits from Theorem 2.3.8 (I). Due to the stochastic self- 
similarity the randomness of the second limit appears only in form of the random 
variable which does not depend on k. In view of Biggin's theorem 2.2.1, M M 
does not vanish with positive probability if and only if 



E 



( E r?) ln + ( E r? 

V i€7i J V i€7i 



< oo 



and in this case EM^ = 1. Then the expectation of the second limit agrees with the 
first limit, i.e. with the limit of the expectations. The integral expression provides 
a formula for numerical calculations in some special situations. Examples for the 
deterministic case may be found in Winter [25] . 

In view of the classical notions the second limit in Corollary 2.3.9 will be called 
random fractal Lipschitz-Killing curvature of order k of the self-similar random set 
F, if it is not zero. 

If the first limit vanishes, one has to check the correctness of the choice of the 
rescaling exponent D — k in order to keep the curvature interpretation. In 'non- 
fractal' situations the exponent has to be replaced by 0. For a detailed discussion 
of this problem see [25] . (Perhaps the rescaling exponents can be used in order to 
distinguish between 'fractal' and 'non-fractal' self-similar sets.) 



3. Proofs of the main results 



3.1. Proof of Theorem 2.3.8. Recall that we wish to reduce the convergence 
assertions in (I) and (II) to the above renewal theorem. Substituting e := R\ J|e~* 
and r := R\ J\e~ s under the integral we obtain the equivalent limit relationships in 
Theorem 2.2.3 for the above introduced processes 

7 ^*C fe ( J F( J R|J'| e -*))l[o,oo)(*) 



Z t 



{R\J\)-"e- 



ieTi 



with Wi = | In r j | . (The use of the constant R\J\ will be clear later.) Therefore it 
suffices to check the conditions on the process Y. 

The measurability and continuity properties in Condition 2.2.2 for the process Y 
follow from its definition together with Lemmas 2.3.1 and 2.3.4, and Corollary 2.3.5. 
For (21) it is sufficient to find some 5 > such that 

- Dt \Y t \ 

ess sup 



E 



t>o 



-St 



< 00 . 



According to the above substitution this may be reformulated as 



(29) 



E 



\Q(r)\ \ 
ess sup [ , n , , I 

0<r<fl|J| \r*-D+S 



J 



< oo . 
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where 



Q(r) := C k (F(r))-Y,r!l ( o,R\J\]£)C k (F i (I r ) 
ieTi n V Tl 

= C k (F(r)) - l(o,ii|7|n]W^(^W) == Mr) 



The equality between Q(r) (where the order k of the curvature is suppressed in the 
notation) and R k (r) (which stands in the assertion of the theorem) follows from 
the scaling property of C k . (Recall that F, t = F n Ji, 3 L = St J , F i = S^ 1 (F) n J, 
and F(r) C U ieTl Mr)-) 
We decompose Q(r) into 

Q(r) = ^ra-^rfC^FQjj 

+ E^Wi,oo)©^(^(ii)) 

ie7i 1 V r * y 

=: Q 2 (r) +Qi(r). 

Using that the F l , i £ 7i, are independent of 71 and have the same distribution as 
-F we obtain for Qi the estimate 



E 



ess sup 

0<r<R\J\ 



(\Qi(r)\\ 


< Ez/E 


\ r k J 





ess sup . 

r>R\J\ \ 



f\C k (F(r))\ 



According to Theorem 4.1 in the Appendix the last expression is finite. Therefore 
it remains to prove that 



(30) 



E 



, IQ 2 (r)| 
ess sup -= — tt-y 

<r<R\J\ \r k ~ D + S 



< oo 



for some S > 0. By the scaling property of C k we get 

(31) Q 2 (r)=C k (F(r))-^C k (F i (r)). 

ieTi 

Next we decompose the total fcth curvatures by means of the corresponding 
curvature measures: 



C k (F(r)) 



C k (F(r), |J Fi(r)\ = C k (F(r),A r U (A r ) c 

^ ieTi ' 
C k (F(r),A r ) + C k (F(r),(A r ) c ), 

A r := [J Jj(r)r\J k {r). 

&k 



where 
(32) 

Similarly, 

C fc (*i(r)) = C k (Fi(r), A r ) + C fc (F,(r), (A r f) , i e Ti . 
The locality of the curvature measure Cfc implies 

Ck(Fi(r), (A r ) c ) = C k (Fi(r), B i ) = C k (F(r), B i ) 
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and F(r) n (A r ) c is the disjoint union of the sets B l := Fi(r) \ A r , i G 7i- Hence, 
C k (F(r), (A r ) c )) - ]T CkiFiir), (A r ) c ) = . 

Substituting this in (31) we infer 

Q 2 (r) = C fe (F(r),A r ) - ^ C fc (Fi(r),^ r ) 

ieTi 

and by the scaling property of Cfc from this 

Q 2 (r) - ft^W^^-^rfftfF^-),^ 1 ^)) 



= Os(r) - Q 4 (r) - Q 5 (r) 



for 



Qs(r) 
Q 4 (r) 



C fe (F(r),A r ) 

E^ 1 (o^|./|](f)^(^(f)^r 1 (A)) 



ieTi 



Q 6 (r) := E^Wl^ftCfcf^&.Sr 1 ^)). 

Therefore it suffices to prove the estimate (30) for Q 3 , Q4, and Q5 instead of Q2 

separately. 

First we obtain 



E 



ess sup 

0<r<fl|J| 



< 00 . 



Here the arguments are the same as for Q\ above taking into regard that for any 
Borel set B, \C k (F,B)\ < C% ar (F,WL d ) and applying Theorem 4.1 in the Appendix. 
For estimating Q 3 and Q4 we we will use the set inclusions 

A r c(SJ) c (r), ^(^nF^-jcJ^-), and J c (r) C (SJ) c (r) , 

(Recall that SJ = UieTi anc ^ ^ s f rom OSC.) Then for Q 3 the estimate follows 
from (33) in Lemma 3.1.1 below. 

Finally, using once more that the sets F l , i <G 7i, are independent of 71 and have 
the same distribution as F, we infer for the same S as above 



E 



ess sup 

0<r<R\J\ \ r 



\Q*(r)\ \ 

,k-D+5 



< E 



r 



k-D+S 



„k-D+5 °* 



< Ei/ E 



< Ei/ E 



ess sup 

0<r<R\J\ 



ess sup 

0<r<il|J| 



o<^<fl|j| y 

'CT(F(r),Jc(r)) 



Ct ar [ FU-),J C (- 



< CO 



according to (33). 



□ 
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Lemma 3.1.1. Under the conditions of Theorem 2.3.8 we have 

'cr(F(r),(Sjy(r)) 



(33) 



E 



ess sup 

0<r<R\J\ 



r k-D+8 



< oo 



for some < S < D. 

Proof. We first proceed similarly as in the above proof choosing the subtree T(r) 
from (26) instead of 71 in the decomposition of the curvature measures. Recall 
that F(r) = UcreT(r) F a {r) for any r > 0. Therefore the locality of the curvature 
measures implies 
(34) 

Cr(F(r),(Sjy(r)) =Cr(F(r),( \J /^(r))n(SJ) c (r)) < £ C? r (F(r), F„{r)) , 

where the subtree %(r) is defined in (28). For k G {d — l,d} wc can use the 
(in) equality 

(35) Cr(F(r),F a (r)) < C k (F a (r)) , a G T b (r) . 

Furthermore, using that the curvature measures are concentrated on the boundary 
of F(r) we obtain for any a G Tb{r) and k < d — 2, 

Cr(F(r),F a (r)) 



= Cr(F(r),F a (r)\ {J F^r)) + C" k aT (V(r), F CT (r) n (J Mr) 
= Cr(F(r),F„(r)\ (J F a ,(r)\ + Cl ar (f{t\ dF a {r) n d( (J F CT ,(r)) 



< 



Cr(^(r))+Cr(F(r),^(r)nS( |J F„,(r))) 

^ a'eT(r), cr'^cr ' 

f k „Cr(F°(f )) +Cr(F(r),dF„(r)nd( \J F a ,(r))). 



Recall that the sets F a — S~ 1 (F) n J, <r G %(r), are independent of the subtree 
7b (r) and have the same distribution as F. Moreover, for a G % (r) we have 
> R\J\- Therefore the above estimates yield 



E 



esssup sup (r- k C% ar (F(r),F a (r))) 

0<r<R\J\ aeT b (r) 



< E 



esssup (r- k C v k ar (F(r))) 

r>R\J\ 



+E 



ess 
o< 



;ssup sup \r- k Cr r (F(r),dF a (r)nd( M F„,(r))) 

r<R\J\ae%(r)\ \ a , eT[r) ^,^ 



The first summand on the right side is bounded by Theorem 4.1 in the Appendix, 
since | J| > (This holds also for k G {d — l,d}.) The boundedness of the second 
summand is assumption (iv) in Theorem 2.3.8. In view of this, (34) and (35) it 
suffices now to show that 



(36) 



E 



sup r D - s $(%(r)) 

0<r<R\J\ 



< CO 
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for some < S < D, where (j denotes the number of elements of a finite set. 

For we use some ideas from the deterministic case (cf. [25] ) and show a probabilistic 

version for our random subtrees 7b(r): 

By the strong open set condition P(Fn int(J) ^ 0) > there exist some a > and 
< p < R\ J\ such that for the subtree 

(37) T{p, a) := {a e T(p) : d{x, dJ)>a, xe F a } 
we have 

(38) P{T(p, a) ^ 0) > . 

The Markov stopping property of the subtree T{p) of T implies 

4 E f?" 

L ^er(p) • 

for the Hausdorff dimension D. Hence, there is a unique < 5 < D such that 



(39) 



E 



E »■■ 

<7€T(p)\r(p,a) 



D-5 



1 . 



For r > denote 

5(r) := {cr = (Ti . . .cr n G T(r) : a k+1 a k+2 T ai - ak (p, a) for any 1 < k < I < n} 

(Recall the notation X a {uj) = X{6 a ijj) for a random clement X.) 

As in the deterministic case (see [25], proof of Lemma 5.4.1, Part I, where the 

polyconvex setting is not needed) one shows that 



iKT 6 (r))<E«(sV)) 



)=i 



where r* := 2{ar m inR\ J\) 1 r. Consequently, 



E 



sup 

0<r<R\J\ 



sup 

r>0 



r d - S m{r*)) 



since the S 1 have the same distribution as 5 and are independent of v. 
Thus it is sufficient to prove that 



(40) 



E 



sup 

r>0 



tt(S(r)) 



< oo . 



Because of the open set condition on the (deterministic) basic set J for r > p' > 
the number of its smaller copies under the random similarities of size of order r 
and hence, the number of elements of T(r) , is uniformly bounded. (Use a volume 
comparing argument for disjoint open balls, one in each copy, of the same radius.) 
Since jj(S(r)) < tt (T(r)) we obtain for any < p' < p a constant C > such that 



sup 

r>p' 



r D - 5 tt(~(r))l <Cw.p.l 



and from this 



(41) 



sup 

r>0 



|(S(r)) ) < max ( sup 

0<r<p' 



«(S(r)) ,C w.p.l. 
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On the other hand, for < r < p the definition of the random set S(r) implies 

m(r))< E «( sff (f))- 

Combining this with (41) for p' := p we infer for any < r' < p, 



*(r') := E 

= E 

< E 

< E 

< E 



max ( sup ( r D ~ 5 Jt(H(r)) ),C 

r' <r 



max( sup [r D - 5 i(E(r))),C 

r'<r<p 



E 



„D-5 



max | sup . , , a _ D _ S 

r '< r <r\*eT( P )\T( P , a) ' 



E 

y€r(p)\T(/p,a) 

E 

Ter(rt\r(p,a) 



rf* 5 max 



s <?<,(^< s ' ( ^)> c 



< E 



E 



f?- 5 max( sup (r D - 5 fl(S»)),C 

*eT(p)\T(p,a) 



r'6<r 



max [ sup ( r 

r'S<r 



tt(H(r)) ,C 



= #(r'0) 



with 9 := p 1 i?|J| > 1, since for a e T(p) we have R\J\r a < p, the random 
numbers tt(S CT (r)) , < r < p, are independently of T(p) \ T(p,a) distributed as 
K(3(r)), and 

_<rer(p)\T(p,a) 

according to (39). By monotonicity of the function vE'(r') we infer for any < r' < p, 
\[r(r') = *(0r') and hence, lim r ^ = which leads together with (41) to 

assertion (40). □ 



4. Appendix - An estimate for Lipschitz-Killing curvature measures 

of large parallel sets 

Here we will show that the variation of the k th Lipschitz-Killing curvature of 
the parallel set of amount r of an arbitrary compact set K for sufficiently large r is 
bounded from above by a constant multiple of r k : 



Theorem 4.1. For any R > \[2 and k = 0, 1,. . .,d there exists a constant Ck(R) 
such that for any compact set K C M d we have for any r > R\K\, 

re&ch(K(f)) > \K\VR 2 - 1, 
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dK(r) is a (d — 1)- dimensional Lipschitz submanifold, and 

C v k ar (K(r),R d ) 

sup k }' '- < c k (R) . 

t>r\k\ r k 

(It is well-known that for compact convex sets K these properties hold for all 
R > 0, the last one in the sharper version of an isodiametric inequality with an 
optimal constant.) 

Proof. According to an argument of Rataj (cf. Lemma 2.2 in Hug, Last and Weil 
[10]) we have for r > \K\, 



reach (K(r)) > yjr 2 - \K\ 2 , 
which implies the first assertion. 

By the scaling property of the curvature measures we get for any r > 0, 

C v k ar (K(r),R d ) = r k C v k ar {{r- 1 K){1), R d ) . 
Hence, it suffices to show for R > y/2 the inequality 
(42) sup C v k ar (K(l),R d ) <c k (R) 

\K\KR- 1 

and that for such K the parallel set K(l) has a Lipschitz boundary. For, we fix an 
arbitrary < s < vT — R~ 2 — R . Since the curvature measures are translation 
invariant, we may assume that € K . Then we obtain for closed balls B(0, r) with 
center and radius r the set inclusions 

5(0, 1 - s) C K(l) C B{0, 1 + \K\) C B{0, 1 + R' 1 ) 

and from this the following estimates for the Euclidean distance function d(-, ■). 

s< inf d(x,K(l)) . 

xedB(0A-s) v ' 

(If d(x,K(X)) < s for some x G 95(0,1 - s), we get d(0,K{l)) < d(0,x) + 

d(x, K(X)} < 1 — s + s = l which is a contradiction, since G K.) Furthermore, 

d(x, KAY)) < d(x, dB(0, 1 + R- 1 )) < 1 + R- 1 - (1 - s) = R- 1 + s, for any x G 
dB(0, 1 — s), which implies 

sup d(x, kJV)) < R- 1 +s< y/l - R- 2 < reach {K{Vj) . 

xedB(0,l-s) 

Therefore we can define a biunique mapping from the unit normal bundle of the ball 
B(0, 1 — s) onto that of the complementary set K{1) as follows: Jk ■ nor£>(0, 1 — 
s) — > nor_ft'(l) with 

, / x\ ( x - u m x 



m) '\x-n W) x\ 



fx is a Lipschitz mapping whose Lipschitz constant is uniformly bounded in K 
as above by some constant c(i?, s) depending only on R and s. This follows from 
Theorem 4.8 in [3]. (Federer's value r is our R^ 1 + s, his q is our — R~ 2 , and s 
has the same meaning as there.) In particular, K(l) has a Lipschitz boundary. 
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The transformation formula for rectifiable currents with compact support under 
Lipschitz mappings (cf. Federer [4, 4.1.30]) implies for the unit normal cycles 

N m (<p)= J <P= J (fK)*<P = N B{0il _ a) ((/x)V) 

norKXT) norB(0,l-s) 

for any smooth differential (d — l)-form ip on R d x M d , where (fx)^ denotes its 
pullback form under the mapping fx (in the sense of a.e. differentiation) . 
Recall that for any set X of positive reach and any bounded Borel set B C K d we 
have 

Ck(X,B) = N X Ll BxR d((p k ) 

for the fcth Lipschitz-Killing curvature form fk, k = 0, . . . , d — 1. This leads to the 
estimates 

C v k ar {K (l),R d ) = Cl ar (K(X),M d ) 

< (Lipt/^^IKorB^i-s)!! \\^(B(0, 1 + R- 1 ) x R d )\\ 

< ciR^y-'WN^s^^W |Kl(B(0, 1 + R- 1 ) x R d )|| 

for the mass norm of rectifiable currents and the comass norm of differential forms 
restricted to B(0, 1 + i? _1 ) x M. d . (For notations and details on current theory see 
Federer [4, Chapter 4].) The minimum over s as above on the right hand side (but 
also the expression for fixed s) provides a desired upper bound Ck(R) if k < d — 1. 
The case k = d is trivial. □ 
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